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Abstract 

^ , We consider a polyharmonic operator H = (—A)' + V(x) in dimension two with 

\q [ I > 2, I being an integer, and a quasi-periodic potential V(x). We prove that 

the spectrum of H contains a semiaxis and there is a family of generalized eigen- 

p . functions at every point of this semiaxis with the following properties. First, the 

eigenfunctions are close to plane waves e^ k ' x ^ at the high energy region. Second, the 
isoenergetic curves in the space of momenta k corresponding to these eigenfunctions 
have a form of slightly distorted circles with holes (Cantor type structure). A new 
method of multiscale analysis in the momentum space is developed to prove these 

!— 1 1 results. 

>- , We study an operator 

H = (-A)' + V(x) (1) 

in dimension two, where I is an integer, I > 2, V(x) is a quasi-periodic potential being a 
trigonometric polynomial: 

O ■ 

(N ; V = V SuS2 e 2ni{si+aS2 ' x) , 1 < Q < oo. (2) 

T"! 1 si,s 2 gZ 2 ,0<|si| + |s 2 |<Q 



We assume that the irrationality measure /x of a is finite: /i < oo, or in other words, that 
a is not a Liouville number 0. 

The one-dimensional situation d — 1, I — 1 is thoroughly investigated in discrete 
and continuum settings, see e.g. PQ-P3] and references there. It is known that a 
one-dimensional quasi-periodic Schrodinger operator demonstrates spectral and trans- 
port properties which are not close to those of a periodic operator. The spectrum of the 
quasi-periodic operator is, as a rule, a Cantor set, while in the periodic case, it has a band 
structure. In the periodic case the spectrum is absolutely continuous, while in the quasi- 
periodic case, it can have any nature: absolutely continuous, singular continuous and pure 
point. The transition between different types of spectrum can happen even with a small 



1 Note, that ji > 2 for any irrational number a. 
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change of a coefficient in a quasi-periodic operator [8]. The mechanism of the difference 
in spectral behavior between periodic and quasi-periodic cases can be explained by a phe- 
nomenon which is known as resonance tunneling in quantum mechanics. It is associated 
with small denominators appearing in formal series of perturbation theory. Since the 
spectrum of the one- dimensional Laplacian is thin (multiplicity 2), resonance tunneling 
can produce an effect strong enough to destroy the spectrum. If a potential is periodic, 
then resonance tunneling produces gaps in the spectrum near the points A n = (irn/a) 2 , 
n G Z, a being the period of the potential. If the potential is quasi-periodic, then it 
can be thought as a sort of combination of infinite number of periodic potentials, each of 
them producing gaps near its own A n -s. Since the set of all A„-s can be dense, the number 
of points surrounded by gaps can be dense too. Thus, the spectrum gets a Cantor like 
structure. The properties of the operator in the high energy region for the continuum case 
d=l are studied in [7J- The KAM method is used to prove absolute continuity of 

the spectrum and existence of quasiperiodic solutions at high energies. 

There are important results on the density of states, spectrum, localization concerning 
the quasi-periodic operators in Z d and, partially, in M. d , d > 1, e.g. [T4]-|24]. However, it 
is still much less known about (JT|) then about its one- dimensional analog. The properties 
of the spectrum in the high energy region, existence of extended states and quantum 
transport are still wide open problems in the multidimensional case. 

Here we study properties of the spectrum and eigenfunctions of ([T]) in the high energy 
region. We prove the following results for the case d = 2, I > 2. 

1. The spectrum of the operator ([1]) contains a semiaxis. 

This is a generalization of a renown Bethe-Sommerfeld conjecture, which states 
that in the case of a periodic potential, / = 1 and d > 2, the spectrum of ([I]) 
contains a semiaxis. There is a variety of proofs for the periodic case, the earliest 
one is [25]. For a limit-periodic periodic potential, being periodic in one direction, 
the conjecture is proved in [26]. For a general case of limit-periodic potential the 
conjecture is proven in [27]- [29]. Here we present the first proof of (a generalized) 
Bethe-Sommerfeld conjecture for a quasi-periodic potential. 

2. There are generalized eigenfunctions \I/oo(k, x), corresponding to the semi-axis, which 
are close to plane waves: for every k in an extensive subset Soo of IR 2 , there is a 
solution ^oc^k, x) of the equation H^/^ = Aqo^/qo which can be described by the 
formula: 

vl/ 00 (k,x) = e^(l+ Moo (k,x)), (3) 
IWUcO*') , = 0{\k\-^), 7i >0, (4) 

k|— i>oo 

where Woo(k, x) is a quasi-periodic function, namely a point-wise convergent series 
of exponentials e J ( n + am ' x > ) n, m G Z 2 . The eigenvalue Aoo(k), corresponding to 
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\l/oo(k, x), is close to \k\ 2L . 

A 00 (k) = |k| 2i + 0(|k|^), 72 >0. (5) 

|k|— 5>oo 

The "non-resonant" set Soo of vectors k, for which (EJ) - (J5]) hold, is an extensive 
Cantor type set: = n™ =1 3 n , where {Snj^Li i s a decreasing sequence of sets in 
M 2 . Each S„ has a finite number of holes in each bounded region. More and more 
holes appear when n increases, however holes added at each step are of smaller and 
smaller size. The set 9oo satisfies the estimate: 

ISooHBrI = |b r |(i + o(ir*)), 73 >o, (6) 

where B R is the disk of radius R centered at the origin, | • | is the Lebesgue measure 
in R 2 . 

3. The set D^X), defined as a level (isoenergetic) set for A QO (k), 

2)oo(A) = {k G Soo : Aoo(k) = A}, 

is proven to be a slightly distorted circle with infinite number of holes. It can be 
described by the formula: 

©00(A) = {k : k = x 00 (A, V)v, V e < B 00 (X)} , (7) 

where B^A) is a subset of the unit circle Si. The set B^A) can be interpreted 
as the set of possible directions of propagation for almost plane waves ([3]). The set 
(A) has a Cantor type structure and an asymptotically full measure on Si as 
A —7- 00: 

L(Boo(A)) = 2tt + 0(A-^), 74 >0, (8) 

here and below L(-) is a length of a curve. The value x oc (A, z7) in ([7j) is the "radius" 
of Dq^A) in a direction V. The function x^A, z7) — X 1 / 21 describes the deviation of 
Doo (A) from the perfect circle of the radius X 1 / 21 . It is proven that the deviation is 
asymptotically small: 

*bo(A,i7) = A 1/2 ' + 0(A^ 5 ), 75 >0. (9) 

A— >oo 

4. The branch of the spectrum corresponding to ^/^(k, x) (the semiaxis) is absolutely 
continuous. 

To prove the results listed above we suggest a method which can be described as 
multiscale analysis in the space of momenta. This is a development of the method, which 
is used in [27]-[29] for the case of limit-periodic potentials. The essential difference is that 



3 



in [27] [29] we constructed a modification of KAM method, where the space variable x 
still plays some role (e.g. in the uniform in x approximation of a limit-periodic potential 
by periodic ones), while in the present situation all considerations are happening in the 
space of the dual variable k. The KAM method in [2 7] -[29] is motivated by |30j-|32j. 
where the method is used for periodic problems. Multiscale analisys which we apply here 
is deeply analogous to the original multiscale method developed in [33] (see also [19], [20] ) 
for the proof of localization. The essential difference is that in [33], [19], [20] the multiscale 
procedure is constructed with respect to space variable x to prove localization, while we 
construct a multiscale procedure in the space of momenta k to prove derealization. 

Here is a brief description of the iteration procedure which leads to the results described 
above. Indeed, let k G K 2 . We consider a set of finite linear combinations of plane waves 
e i(k+p+am,x)^ p, m G Z 2 . The set is invariant under action of the differential expression 
(PQ). Let HQs.) be a matrix describing action of ([I]) in the linear set of the exponentials. 
Obviously, 

H(k) = H (k) + V, #o(k)(p,m),(p',m') = |k + p + am|f 2 (5(p,p')(5( m>m /), 

^(p,m),(p',m') ^p— p',m— m'- 

Next, we consider an expanding sequence of finite sets M n in the space Z 2 x Z 2 of indices 
(p,m): M n C M n+ i, lim^oo M n = Z 2 x Z 2 . Let P n be the characteristic projection of 
set M n in the space £ 2 (Z 2 x Z 2 ). We consider a sequence of finite matrices H^ n \k) = 
P n H(k)P n . Each matrix corresponds to a finite dimensional operator in £ 2 (Z 2 x Z 2 ), given 
that the operator acts as zero on (I — P n )£ 2 . For each n we construct a "non-resonant" set 
S„ in the space M 2 of momenta k, such that: if k G Sn> then H( n >("k.) = P n H(k)P n has an 
eigenvalue A n (k) and its spectral projector £ n (k) which can be described by perturbation 
formulas with respect to the previous operator i7^ n_1 ^(k). If k G nJ^Sn then A n (k) 
and £ n (k) have limits. The linear combinations of the exponentials, corresponding to the 
projectors £ n (k), have a point-wise limit in x, the limit being a generalized eigenfunction 
of ([!]) . The generalized eigenfunction is close to the plane wave e* <k ' x> in the high energy 
region. 

Each matrix is considered as a perturbation of a matrix H^ n \ the latter has a 
block structure, i.e., consists of a variety of blocks (k + p + am) , s = 1, ...,n — 1, and, 
naturally, some diagonal terms. Blocks with different indices (s) have sizes of different 
orders of magnitude (the size increasing with s). Thus we have a multiscale structure in 
the definition of H^ 1 '. We use H^(k) as a starting operator to construct perturbation 
series for H^(k). At a step n we apply our knowledge of spectral properties of H ^ (k + 
p' + ara'), s = 1, ...,n — 1, p', m' G Z 2 , obtained in the previous steps, to describe spectral 
properties of H^(k + p + am), p, m G Z 2 and to construct 9n- 

At step one we use a regular perturbation theory and elementary geometric consider- 
ations to prove the following results. There is a set 9i C M 2 such that: if k G Si, then 
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the operator H^(k) has a single eigenvalue close to the unperturbed one: 

A«(k) = |k| 2i + 0(|k|-^), 72 >0. (10) 

k|— >oo 

A normalized eigenvector u« is also close to the unperturbed one: u*- 1 -* = u*- -* + u^, 
where (u^)( Piin ) = 5 p0 5 m0 and the / 1 -norm of uW is small: ||uW||/i < |k|~ 71 , 7l > 0. It 
follows that: 

v& 1 (k,x) = e^+n 1 (k ) x), \\ui\\l^) , = 0(|k|^), Tl > 0, (11) 

|k|— >oo 

where ^(k, x), Wi(k, x) are the linear combinations of the exponentials corresponding to 
vectors and u^, respectively. It is shown that function \l/i(k, x) satisfies the equation 
for eigenf unctions with a good accuracy: 

-A*i + V* 1 = \k\ 2l * 1 + f 1 , H/ilkoCR*) , = 0(|k|-^) 76 >0. (12) 

|k|— >oo 

Relation (Tl0|) is different iable: 

VA (1) (k) = 2/|k| 2 '- 2 k + 0(|kr 77 ), 77 >0. (13) 

k|— >oo 

Next, we construct a sequence Sn, n > 2, such for any k G S« the operator H^(k) has a 
single eigenvalue A^(k) in a super exponentially small neig hborhood of A< n-1) (k): 

A (n) (k) = A ( "- 1) (k) + 0(|kr |k|78 "), 78 >0. (14) 

|k|— >oo 

Similar estimates hold for the eigenvectors and the corresponding functions \l/ n (k, x): 
#„(k,x) = *„-i(k,x) + u n (k,x), \\u n \\ Looim = O (|k|~ |k|79n ) , 79 >0. (15) 

|k|— >-oo 

-A*„ + y* n = A^(k)* n + / n , 11/nlk.c*) , = O (|k|-l k l 710 ") , 7l0 >0. (16) 

k|— 5>oo 

Formula (Tl4|) is different iable with respect to k: 

VA (n) (k) = VA (n - 1) (k) + 0(|k|-l k l 78 "), 78 >0. (17) 

|k|— >oo 

In fact, for large n estimates (|14j) - (ITTj) are even stronger. The non-resonant set Sn is 
proven to be extensive in M 2 : 

|3„nB R | = \B K \(1 + 0(R-™)). (18) 
Estimates (fl4|) - (fl~8|) are uniform in n. 
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Figure 1: Isoenergetic curve ©i(A) 



Figure 2: Isoenergetic curve 2) 2 (A) 



The set D n (X) is defined as the level (isoenergetic) set for the non- resonant eigenvalue 
A< n >(k): 

D n (A) = {keS„:AW(k) = A}. 

This set is proven to be a slightly distorted circle with a finite number of holes (see Fig. 
Q], [2]). The set 2)„(A) can be described by the formula: 

D n (A) = {k : k = x„(A, z7)z7, z7 G £ n (A)} , (19) 

where !B n (A) is a subset of the unit circle Si. The set 23 n (A) can be interpreted as the set 
of possible directions of propagation for almost plane waves \I/ n (k, x), see (fTTj) . f|T5|) . It 
has an asymptotically full measure on S'i as A — > oo: 

L(S n (A)) = 2tt + 0(A^ 4/2 '). (20) 

A— ^oo 

Each set S n (A) has only a finite number of holes, however their number is growing with 
n. More and more holes of a smaller and smaller size are added at each step. The value 
x n (X, v) — A 1 / 2 ' gives the deviation of D n (X) from the perfect circle of the radius A 1 / 2 ' in 
the direction v. It is proven that the deviation is asymptotically small: 

x n (A,z7) = A 1 / 2i + 0(A^), dKn{ ^ V) =0(X-^), 75,711 >0, (21) 

<p being an angle variable, V = (cos (p, sin <p). Estimates fl20|) . fl2TT) are uniform in n. 

On each step more and more points are excluded from the non- resonant sets S„, 
thus {Sn}^i is a decreasing sequence of sets. The set Soo is defined as the limit set: 
Soo = n^ x 9n- It has an infinite number of holes, but nevertheless satisfies the relation 
(J6j). For every k G Soo and every n, there is a generalized eigenfunction of H^ n ' of the 
type (fTTl) . ( fTol) . It is proven that the sequence of ^ n (k, x) has a limit in /^(IR 2 ) when 
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k G Soo- The function \I/oo(k, x) = lim^oo \l/ n (k, x) is a generalized eigenfunction of H. 
It can be written in the form (J3j) — (SD - Naturally, the corresponding eigenvalue Aoo(k) is 
the limit of A^(k) as n — > oo. 

It is shown that {!B n (A)}^L 1 is a decreasing sequence of sets, on each step more and 
more directions being excluded. We consider the limit 23 00(A) of 23 n (A): 

00 

IMA) = p| 23 n (A). (22) 

n=l 

This set has a Cantor type structure on the unit circle. It is proven that 23 00(A) has 
an asymptotically full measure on the unit circle (see (JSJ)). We prove that the sequence 
x n (A, z7), n = 1, 2, describing the isoenergetic curves D n (A), quickly converges as n — > 
00. We show that Doo(A) can be described as the limit of D n (A) in the sense (J7|), where 
Xoo(A, v) = limn-j.oo x n (A, v) for every V G 23 00 (A). It is shown that the derivatives of the 
functions x n (A, v) (with respect to the angle variable on the unit circle) have a limit as 
n — > 00 for every z7 G 23 00 (A). We denote this limit by d * ao J^' u *) . Using ( 12"T|) . we prove that 

9xoo(A, z7) 



<9<£> 



O (A -711 ) . (23) 



Thus, the limit curve Doo(A) has a tangent vector in spite of its Cantor type structure, 
the tangent vector being the limit of corresponding tangent vectors for D n (X) asn-> 00. 
The curve 2) 00 (A) looks as a slightly distorted circle with infinite number of holes for every 
sufficiently large A, A > A*(V). It immediately follows that [A*, 00) is in the spectrum of 
H (Bethe-Sommerfeld conjecture). 

The main technical difficulty to overcome is the construction of non-resonant sets 
23 n (A) for every fixed sufficiently large A, A > \q(V), where Ao(V^) is the same for all n. The 
set 23 n (A) is obtained by deleting a "resonant" part from 23„_i(A). Definition of 23 n _x(A) \ 
23 n (A) includes eigenvalues of H^-^ik). To describe 23„_i(A) \ 23„(A) one has to consider 
not only non-resonant eigenvalues of the type (flOl) . ([ED , but also resonant eigenvalues, 
for which no suitable formulas are known. Absence of formulas causes difficulties in 
estimating the size of 23 n _x(A) \ 23 n (A). To treat this problem we start with introducing 
an angle variable ip G [0,27r), V = (cosy?, sin^) G Si and consider sets 23„(A) in terms of 
this variable. Next, we show that the resonant set 23 n _i(A) \ 23 n (A) can be described as 
the set of zeros of functions of the type 

det(# (s) (x n _i(^) + p + aim) -A-e), s = l,...,n-l, (p, m) G M n \ (0, 0), 

where x n _i(<£>) is a vector-function describing D n _i(A): x n _i((^) = x n _i(A, u)u. To obtain 
23„_i(A) \ 23 n (A) we take all values of e in a small interval and (p, m) in some subset of 
M n . Further, we extend our considerations to a complex neighborhood $0 of [0, 2ir). We 
show that the determinants are analytic functions of ip and, by this, reduce the problem of 
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Figure 3: Set <P 2 




Figure 4: Set <P 3 



estimating the size of the resonant set to a problem in complex analysis. We use theorems 
for analytic functions to count zeros of the determinants and to investigate how far the 
zeros move when e changes. It enables us to estimate the size of the zero set of the 
determinants, and, hence, the size of the non-resonant set (P n C $o, which is defined as 
a non-zero set for the determinants. Proving that the non-resonant set (£> n is sufficiently 
large, we obtain estimates (TIB"]) for Sn and ( 1201) for !B n , the set !B n corresponding to the 
real part of <P n . To obtain <P n we delete from $o more and more discs (holes) of smaller 
and smaller radii at each step. Thus, the non-resonant set <P n C $o has a structure of 
Swiss Cheese (Fig. [31 H]). Deleting a resonance set from <P at each step of the recurrent 
procedure we call a "Swiss Cheese Method" . The essential difference of our method from 
constructions of non-resonant sets in similar situations before (see e.g. (3D] [32]. [31]) is 
that we construct a non-resonant set not only in the whole space of a parameter (k e M? 
here), but also on isoenergetic curves D n (X) in the space of the parameter, when A is 
sufficiently large. Estimates for the size of non-resonant sets on a curve require more 
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subtle technical considerations than those sufficient for description of a non-resonant set 
in the whole space of the parameter. But as a reward, such estimates enable us to show 
that every isoenergetic set for A > Ao is not empty and thus, to prove Bethe-Sommerfeld 
conjecture. 

Note that generalization of the results from the case / > 1, I being an integer, to the 
case of rational I satisfying the same inequality is relatively simple; it requires just slightly 
more careful technical considerations. The restriction / > 1 is also technical, though it 
is more difficult to lift. The condition I > 1 is needed only for the second step of the 
recurrent procedure. The authors plan to consider the case Z = 1 in a forthcoming paper. 
The requirement /i < oo is essential, since we use it to estimate the minimal values of 
|p + am | when (p, m) e M n \ (0,0). Such estimates are necessary for controlling small 
denominators in the perturbation series at each step. 

Acknowledgement The authors are very grateful to Prof. Leonid Parnovski for useful 
discussions and to Prof. Young- Ran Lee for allowing us to use figures 1-4 from [27J. 
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